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Abstract. We selected a sample of galaxies, extremely isolated in 3D redshift 
space, based on data from NED and the ongoing ALFALFA HI (21cm) survey. 
A simple selection criterion was employed: having no neighbors closer than 300 
km-s -1 in 3D redshift space (a, S, z). The environments of galaxies, selected 
using this criterion and NEEQ data alone, were analyzed theoretically using a 
constrained simulation of the local Universe, and were found to be an order of 
magnitude less dense than environments around randomly selected galaxies. One 
third of the galaxies selected using NED data alone did not pass the criterion 
when tested with ALFALFA data, implying that the use of unbiased HI data 
significantly improves the quality of the sample. 



Introduction 



The research described in this paper is part of an extensive study of star forma- 
tion and evolution of galaxies in different environments and of va rious morpho- 
logical types, conducted by our group f or the past few decades ( Broschl Il983l ; 



Brosch et al.lll998l : IZitrin fc Broschl 12003 ). Specifically, we are studying galaxies 



in the most extremely underdense regions of the local Universe. These galaxies 
are particularly interesting, since they have evolved with little or no environmen- 
tal interference, and therefore are good test beds for calibrating galaxy evolution 
models. Furthermore, when compared to galaxies in denser regions they illumi- 
nate the overall effects of the environment on the evolution of galaxies. 

The first and possibly most crucial step in this study is selecting a sample 
of extremely underdense region galaxies, referred to here as Isolated Galaxies 
(IGs) . By definition these IGs are very rare, and therefore at moderate redshifts 
the sample is fairly small, including only a few dozen galaxies. Here we describe 
the method used for selecting the sample and the validation process used to 
ensure that these are indeed galaxies in the most extremely underdense regions 
of the local Universe. 

The sample is currently being imaged in optical wide bands and H a to mea- 
sure the current star formation and the evolutionary histories of these galaxies. 
The obser vations are carried out at the Wise Observ atory. The ALFALFA blind 
HI survey (|Giovanelli et al.ll2005l ; iMartin et al.ll2009l ). which covers 7000 deg 2 of 
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the sky to redshift 17000 km-s -1 with low detection limit and excellent position- 
ing accuracy (~0.1'), is a key data source for selecting the sample and analyzing 
i ts HI content. Furthe r observational data are obtained from NED and SDSS 
(jAbazaiian et al.ll2009h . 



2. The Isolation Criterion 

In the last few decades great advances were made in redshift surveys, which now 
map the local Universe in redshift space with great precision. Before these were 
available isolated galaxies had to be identified using projected coordinates alone 
(a, 5), i.e. searching in two-dimensional space (2D). Sophisticated algorithms 
were used to estimate radia l distance, e.g. via th e angular sizes of galaxies, 
such as done for the classical iKarachentseval ( 19731 ) catalog of isolated galaxies. 



Nowadays, when the local Universe is mapped in fine detail, it is possible to 
perform three-dimensional (3D) searches. The advantages of using such strategy 
are simplicity and straightforwardness, not having to assume anything about the 
characteristics of the galaxies (size, magnitude, etc.). 

However, using redshift mapping introduces two difficulties. First is the 
incompleteness of most redshift databases. A galaxy that seems to be isolated 
might have neighbors for which a redshift was never measured. Second is the 
error in radial distance, introduced by peculiar velocities. Using redshift data 
one performs a search in 3D redshift space, the mathematical representation 
of the projected (a, 5) and radial (redshift) coordinates. It should be kept in 
mind that mapping in 3D redshift space (a, 5, z) can differ significantly from the 
true mapping in position space (a, 5, r). For example, although close in position 
space, two galaxies in a cluster might have very different redshifts due to cluster 
velocity dispersion and will therefore seem distant in 3D redshift space. 

In this work we have chosen to use a simple selection criterion. A galaxy 
is considered an IG and is included in the sample if it has no known neighbors 
closer than 300 km-s -1 in 3D redshift space, and its redshift is in the range 
2000 < cz < 7000 kirns" 1 . 

The redshift range is limited to 7000 km-s -1 to have better completeness of 
redshift data around each galaxy. The reason for the lower limit of 2000 km-s -1 
is to keep the sky area that has to be searched around each galaxy, relatively low 
(at 2000 km-s -1 neighbors have to be searched for as far as a projected distance 
of 8.6°). 



3. The Sample 

The search criterion was applied to two sky regions, one in the spring sky: 
7h30m< a <16h30m, 4° < 5 < 16°, 2000 < cz < 7000 km-s -1 , and the other 
in the autumn sky: 22h00m< a <03h00m, 24° < 5 < 28°, 2000 < cz < 7000 
km-s -1 . Both regions were already c overed by AFLAFA. T he "spring" sky region 
is also fully covered by SDSS DR7 (jAbazaiian et al.1 12009^ . 



IGs were searched for in these regions first based only on NED data of Jan- 
uary 2009. In the "spring" sky region 27 IGs were found out of 2826 candidates 
(~1.0%). In the "autumn" sky region 6 IGs were found out of 244 candidates 
(~2%). ALFALFA data were then used to find possible neighbors to the "spring" 
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IGs. Out of 27 IGs, nine (a third) had ALFALFA neighbors. Four galaxies had 
more than one ALFALFA neighbor, and one galaxy had 7 ALFALFA neighbors, 
none of which had redshift data from optical spectra! 

The search was repeated for the "spring" sky region in a different fashion. 
The search criterion was first applied to ALFALFAs data, and then the isolation 
of the resulting galaxies was tested using NED. Five additional IGs were found 
this way (none of which had redshift data in NED). 

Figure 1 compares the number density n around a typical IG from the 
sample (UGC 12123) to the number density around a typical field galaxy (LEDA 
166859) and around M87, located at the center of the Virgo cluster. For each of 
the three galaxies n is shown as a function of r, the radius of a sphere around 
the galaxy for which n was calculated. 

UGC 12123 has two neighbors within 5 Mpc-h" 1 : one 3.06 Mpc-h" 1 away 
(detected by ALFALFA) and the other 4.35 Mpc-h" 1 away. As can be seen in 
Figure 1, the neighborhood density of the IG (UGC 12123) is about one order of 
magnitude lower than that of the typical field galaxy (LEDA 166859) and more 
than two orders of magnitude lower than that of the cluster galaxy (M87). 

4. Validation 

In order to validate that when applying the search criterion to NED data, one 
obtains a sample of galaxies in extremely underdense regions of the Universe, 
an independent dataset was used. This dataset was the Box 160 constrained 
simulation of the local Universe (Hoffman & Gottlober 2008, private commu- 
nication; Gottlober S. & Klypin 2008; Forero- Romero et al. 2009). Box 160 is 
a ACDM (WMAP3) based mock universe that emulates large structures in the 
local Universe (Virgo, Coma, Local Supercluster, etc.) with box side length of 
160 Mpc-h" 1 and 1024 3 particles, 2.54-10 8 M -h" 1 each. 

The Box 160 dataset was used to create mock "observable" data, consist- 
ing of (a, 5) coordinates and observable magnitudes for all galaxies. Redshift 
data were added to the mock "observable" data for 38% of the galaxies with 
observable magnitude m < 18.5 (randomly selected) and non of the galaxies 
with m > 18.5. This was done to simulate the dilution of redshift data in the 
NED dataset used for selecting the sample from the "spring" sky region. 

A "mock sample" was created by applying the search criterion to the "mock 
observable" data. The "true neighborhood" of the "mock sample" IGs was 
analyzed using the full Box 160 data, i.e. using "true" coordinates of all the 
neighbors including those which the "mock observable" dataset missed. 

The results are shown in Figure 2, where the likelihood function of the num- 
ber density around a galaxy is plotted for the "mock sample" and for all the Box 
160 galaxies ("all"). As can be seen, the IGs selected by this method (based on 
data that simulate NED) are located in significantly underdense neighborhoods 
compared to the general population. 

As we have seen, the use of ALFALFA significantly improves the quality 
of the sample, i.e. the likelihood function is expected to be shifted to lower 
densities. However, this could not be simulated because a mock universe dataset 
that includes HI mass data is not available at present. Without HI data it cannot 
be discerned which simulated objects would have been observable by ALFALFA. 
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Figure 1. The number density n in a sphere of varying radius r around 
three galaxies: UGC 12123 an IG from the sample, LEDA 166859 a field 
galaxy, M87 the central galaxy of the Virgo cluster. 

Figure 2. Likelihood function of the number density n in a 5 Mpc-h -1 
radius sphere around galaxies of the "mock sample" , and around galaxies of 
the entire mock universe (All). Errors are indicated by the + signs in the 
legend. 

5. Conclusions 

• Redshift surveys provide enough information for choosing IG samples based 
on 3D mapping. 

• The simple isolation criterion (having no known neighbors within 3 Mpc-h -1 ), 
when applied to redshifts in the range 2000 < cz < 7000 km-s -1 , provides 

a sample of galaxies in the low density regions of the local Universe. 

• ALFALFA data are extremely useful in eliminating false positives from 
the sample due to the redshift data it provides for low-luminosity galaxies 
with high HI masses. 
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